Laser-heater assisted plasma channel formation in capillary discharge waveguides 



O 
(N 



N. A. Bobrova,i'2 P. V. Sasorov,^'^ C. Benedetti,^ S. S. Bulanov,^ C. 

G. R. Geddes,"* C. B. Schroeder,'* E. Esarey,^ and W. P. Leemans^-^ 

^Institute of Theoretical and Experimental Physics, Moscow 117218, Russia 
^Moscow Institute of Physics and Technology, Dolgoprudny, Russia 
Keldysh Institute of Applied Mathematics, Moscow, Russia 
'^Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA 
'^University of California, Berkeley, California 94720, USA 

A method of creating plasma channels with controllable depth and transverse profile for the 
guiding of short, high power laser pulses for efficient electron acceleration is proposed. The plasma 
channel produced by the hydrogen-filled capillary discharge waveguide is modified by a ns-scale laser 
pulse, which heats the electrons near the capillary axis. This interaction creates a deeper plasma 
channel within the capillary discharge that evolves on a ns-time scale, allowing laser beams with 
smaller spot sizes than would otherwise be possible in the unmodified capillary discharge. 
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Extending the interaction length of an intense laser prop- 
agating in an underdense plasma by using a preformed 
channel is important for a wide variety of applications, 
including laser-plasma acceleration [l| , x-ray lasers 0, Q i 
harmonic generation and direct laser acceleration 
Efficient laser-plasma accelerator (LPA) operation 
requires laser pulse propagation over a distance of the or- 
der of the acceleration length without significant diffrac- 
tion. Plasma channels, i.e., a plasma column with a min- 
imum plasma density along the laser propagation axis, 
may be employed to optically guide short-pulse lasers, 
thereby extending the interaction length Plasma 
channels may be created using a variety of techni ques , 
including laser-induced hydrodynamic expansion |7Hll|. 
ablative capillaries |12 - [l4| . fast Z-pinch assisted channel 
formation |15| . and gas-filled discharge capillaries (l6l - 
[Tsl . Hydrogen gas-filled discharge capillaries have been 
successfully used to generate electron beams in the GeV 
energy range [l9M2l| . 

Further increase in accelerated electron bunch energy 
lies in the utilization of more powerful lasers interact- 
ing with lower density plasmas. The BELLA (Berkeley 
Lab Laser Accelerator) project will employ a PW-class 
laser and up to a meter-scale plasma channel with on- 
axis density of about 10^^ cm~^ [l^, allowing 10 GeV 
electron bunch energy to be reached, and demonstrating 
the feasibility of laser plasma acceleration as a possible 
technology for the next generation of lepton colliders . 
However while the requirements of the maximal electron 
energy and available laser power fix the density and the 
spot size these can not be satisfied simultaneously by 
the available capillaries. Either the matched spot size is 
larger than required for laser guiding, or the small radius 
capillary, which would support guiding with the required 
matched spot, will experience extensive wall damage pre- 
venting the repetitive utilization. 

In what follows we propose a robust method to guide 
short duration, high power laser pulses with the required 
spot size in a low density plasma channel, generated in 



a hydrogen-filled capillary discharge waveguide. After a 
parabolic density profile is established in the capillary by 
the discharge, a long (ns) laser pulse is fired along the 
capillary axis. The laser heats the electrons near the axis 
creating a secondary channel atop the already existing 
channel created by the discharge. The secondary channel 
width and depth can be controlled by varying the heat- 
ing pulse duration and intensity, thus providing means 
of engineering a channel with the properties required for 
efficient acceleration of electrons. 

In order to study the spatial and time evolution of 
the plasma electron density and temperature in a cap- 
illary discharge we performed computer simulations us- 
ing a one-dimensional two temperature (ion and electron) 
single fiuid magnetohydrodynamic (MHD) code The 
fact that the length of the capillary is much larger than its 
diameter justifies the utilization of the one-dimensional 
approximation, which is backed by the experimental re- 
sults [l^l showing the distribution of the discharge pa- 
rameters along the capillary axis to be uniform. 

The physical model behind the MHD code was used to 
simulate the dynamics of dense Z-pinches (25j and both 
discharge-ablated HI] and pre-filled by gas capillary dis- 
charges [23|. We simulate the alumina capillary with an 
inner radius, Rq, of 500 fim, pre-filled with pure hydrogen 
of uniform temperature and density of 5.0 x 10~^ g/cm^. 
The discharge is initiated by a pulse of current, which is 
approximately sinusoidal, with a half-period of ^ 800 ns 
and a peak value of 400 A, driven by an external circuit. 
The parameters of the capillary and the discharge are 
typical for the experiments with high power laser guid- 
ing. The plasma behavior is shown in Fig. 1. After 
the breakdown the current pulse heats the plasma and 
creates an azimuthal component of the magnetic field. 
However, pinching of the plasma is negligibly weak since 
the plasma pressure is much higher than the magnetic 
field pressure. The plasma is confined from expansion in 
the radial direction by the capillary walls. We see three 
stages of plasma evolution. In the first stage (< 100 ns). 
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the radial distributions of the electron density, ion and 
electron temperatures and the degree of ionization re- 
main homogeneous until the moment of almost full ion- 
ization of hydrogen. The electric current penetrates the 
plasma very quickly, and plasma is heated and ionized 
locally. The redistribution of plasma across the capillary 
channel begins at w 100 ns and lasts until w 250 ns. Dur- 
ing this phase thermal conduction becomes significant. In 
the third stage the plasma remains almost immobile in 
the radial direction. The ion and electron temperatures 
are equal. The radial distribution of plasma parameters 
corresponds to a quasi steady-state equilibrium at a given 
electric current. The plasma pressure is almost homoge- 
neous. The Ohmic heating is balanced mainly by ther- 
mal conduction to the relatively cold capillary wall. The 
temperature has its maximum on the axis and the hydro- 
gen plasma is fully-ionized. Then, at constant pressure, 
a monotonic radially-decreasing temperature results in 
an axial minimum in the electron density profile of the 
plasma. During the interval 250 ns < i < 800 ns the ra- 
dial distribution of electron density is almost constant - 
at least to the extent that longitudinal plasma flow may 
be neglected. The time scale of the outflow is of the or- 
der of L/cs, where L is the length of the capillary, and 
Cs = \jTeJr(ii is the ion acoustic velocity. For L = Im 
the time of the longitudinal plasma outflow will be of the 
order of 30/is. 
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FIG. 1: Capillary discharge dynamics without laser heater: 
The distribution in tlie (r, plane of electron density in units 
oflO^'^cm"" 
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For t > 250 ns there is no screening of the axial electric 
field, and consequently it is uniform across the diameter 
of the capillary. There is no difference between electron 
and ion temperatures. This allows us to use a simple 
model of the plasma under such equilibrium conditions, 
formulated in [23|. According to this model, a Gaussian 
beam will propagate through a plasma channel with a 
parabolic electron density profile of the form ne(r) — 



Woifim] « 1.48 X 10^ 



(zn,o[cm-3])i/4 



(1) 



where Rq is the capillary radius, z is the mean ion charge, 
and UiQ is the initial ion density. For the case under con- 
sideration we find that Wq « 145 /xm, which is typically 
too large for efficient LPA operation in the linear regime, 
i.e., to avoid strong laser evolution at high intensities 
{e.g., 60 — 90 fim is required for BELLA experiment de- 
pending on the regime of operation [2^). In order to 
decrease the matched spot size it is necessary to use a 
capillary of a smaller radius. However, operation with 
small radii results in to high wall lattice temperature 
and to melting or disruption of the wall material both 
by the discharge and the wings of the high power laser 
pulse which is to be guided through this capillary. 
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FIG. 2: Capillary discharge dynamics with laser heater: 
distribution of the electron density in units of 10^^ cm 
the (r, t) plane. 
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We propose to modify the plasma density profile in 
such a way that the capillary will gain the capability to 
guide pulses with smaller spot sizes, i.e., to provide ad- 
ditional heating of the capillary plasma near the axis by 
the absorption of the laser radiation. The heating laser 
pulse is fired along the capillary axis and is implemented 
into the MHD code in the following form: 



P(r,i) = Pocos^ 



/ TTr 

V 2P/as 



At 



(2) 



where P{r, t) is the heating laser pulse power and is equal 
to outside the intervals: ti < t < tf and r < Rias- 
The constant coefficient Pq is determined by the en- 
ergy of laser pulse, £, the parameter Rias and duration 
At = tf — ti. The FWHM duration of the pulse is equal 
to At/2, and Rias = 7rWo/2-\/2, which is chosen to match 
the form ^ of the laser pulse to a Gaussian pulse with a 
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spot size of Wq . We set the following parameters for the 
simulation: £ = 13, At/2 = 1 ns and Rias = 161 fJ-Ta. A 
one micron laser is assumed with critical electron density 
ricr = 1.1 • 10^^ cm~'^. The time of the maximum of elec- 
tric current is chosen as the start of laser heating (t=400 
ns) . At that time the plasma is in equilibrium, hydrogen 
is fully ionized and there is a well established parabolic 
electron density profile to guide the heating laser pulse. 
The duration of the heating laser pulse (~ Ins) and the 
time of its propagation (~ 3ns) in the capillary with 
length L = Im are much less than the hydrodynamic time 
(~ 15ns) at i « 400ns. Due to this separation of time 
scales, the energy deposition of the laser radiation can be 
considered as homogeneous along the capillary. So we can 
consider a ID (in radial dimension) approximation. We 
do not take into account the change of the matched spot 
size during the propagation of the heating laser pulse. 
The evolution of the radial distribution of electron den- 
sity under the action of the heating laser pulse is plotted 
in Fig. 2. During the interval 401 ns < t < 415 ns 
the radial distribution of electron density near the axis 
is different from the equilibrium one. An approximately 
parabolic electron density profile with smaller matched 
spot size exists over several nanoseconds. 

We use the results of Ref. [l^l to calculate the matched 
spot size of the high power laser pulse, guided by a modi- 
fied plasma density profile. In Fig. [3^ the dependence of 
the matched spot size on the delay between the heating 
pulse and the high power laser pulse is shown. One can 
see that the matched spot size oscillates as the delay is 
increased and there are two distinct minima at t ^ 402 
ns and at t ~ 417 ns. These oscillations are due to the 
fact that the cavity generated by the heating laser pulse 
has a form of a bubble (see Fig. 2), where the depth and 
radius change simultaneously as we move further away 
from the heating laser pulse along the bubble central axis. 
Since the matched spot size is the function of the chan- 
nel depth, there are two points along the central axis that 
correspond to the matching of the laser pulse to the mod- 
ified plasma density profile. In Fig. 3a these two points 
correspond to the local minima in i? = R{t). We use 
the first minimum to guide the high intensity laser pulse, 
since it supports guiding of smaller matched spot size 
laser pulses. Although R changes as a function of time 
near the minimum matched spot size, on a nanosecond 
scale the variation is sufficiently small to significantly af- 
fect the guiding. The ns level of synchronization between 
lon g an d short pulses was demonstrated in experiments 

The propagation of the laser pulse with the anticipated 
parameters of BELLA laser (40 J, 100 fs and Wq = 64 
/im) was modeled with the code INF&RNO (INtegrated 
Fluid & paRticle simulatioN cOde), which is a 2D cylin- 
drical code that solves the full wave equation for the laser 
envelope evolution coupled to the equation of motion for 
the plasma (PIC and cold fluid descriptions are avail- 
able) and to Maxwell equations for the electromagnetic 
field in the plasma. The plasma responds to the laser 
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FIG. 3: a) The dependence of the laser pulse matched spot 
size, Wo, on the delay between the guided pulse and the heat- 
ing pulse for different values of ao = 0, 0.5, and 1, with 
(black curves) and without (red curves) laser heating, b) 
The dependence of ao on the propagation distance for the 
high power laser delay of ns (black), 2.5 ns (red), and 23.5 
ns (blue). The dashed line corresponds to a ID simulations. 
Here ao = 1.7. The simulations for blue and black curves 
were stopped at z ~ 23 cm, because in these two cases laser 
pulses show large mismatch oscillations. 



field via the averaged laser ponderomotive force. The 
adoption of the cylindrical geometry allows the descrip- 
tion of 3D physics (laser evolution, electromagnetic field 
structure) at 2D computational cost. The code features 
an improved laser envelope solver which enables an ac- 
curate description of the laser pulse evolution deep into 
depletion even at a reasonably low resolution. Additional 
details on INF&RNO may be found in Ref. H. 

The results are shown in Fig. ^jp, where the depen- 
dence of ao, the peak value of the dimensionless vector 
potential of the laser electromagnetic field, on the prop- 
agation distance for the plasma density profiles at 400 
ns, 402.5 ns, and 423.5 ns is presented. The results of 
ID simulations (the dashed curve) demonstrate that the 
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slow increase of gq is due to the longitudinal pulse evo- 
lution and not to the transverse dynamics. From the 
close matching of the ID uq evolution to the 2D result in 
red in Fig. 3b one can see that modified plasma density 
profile provides matched guiding of a high power laser 
pulse with a 1.5 ns delay after the peak of the heating 
pulse. This delay corresponds to the minimum matched 
spot size that can be achieved with this modified density 
profile. The depth and width of the modified density pro- 
file can be controlled by the energy and duration of the 
heating laser pulse, making it possible the matching of 
different high power laser spot sizes for efficient guiding. 

We notice that the dependence of the peak value of qq 
on the propagation distance demonstrates oscillatory be- 
havior. The reason for the oscillations is connected with 
the fact that it is not possible to ensure the matching 
to the plasma density profile uniformly along the high 
power laser pulse. The oscillations of the non-matched 
parts of the laser pulse as well as their interaction be- 
tween themselves and with the matched part result in 
these non-regular oscillations, seen in Fig. 3b. 

In conclusion, we have proposed a robust method of 
controllable modification of the plasma density profile in 



the gas-filled capillary discharge waveguide. It is achieved 
by heating the plasma near the capillary axis by a long 
(ns scale) laser pulse. The laser heating results in the 
generation of a secondary plasma channel atop of the ex- 
isting discharge-generated channel, which serves the pur- 
pose of guiding high power laser pulses with smaller spot 
sizes than the non-modified plasma density profile can 
provide. We demonstrated through MHD modeling that 
such profiles can be generated and have shown that they 
should be capable of guiding high power lasers. This 
method enables larger capillary radii, thereby prevent- 
ing the extensive wall damage and thus greatly increases 
the lifetime of the capillary. This is extremely impor- 
tant for the LPA operation in the high repetition regime. 
The multidimensional effects (including laser heater evo- 
lution), neglected in the present paper, will be a subject 
of future studies for different capillary geometries. 
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